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ABSTRACT 


: 


Internal  waves  scatter  sound  by  two  related  perturbations:  (1)  those 

associated  with  vertical  particle  displacements  r.(x,y,z,t)  in  the  presence 

of  a vertical  gradient  of  (potential)  sound  speed:  6c  ■ r,*3  c ; (2)  those 

z p 

associated  with  horizontal  particle  velocities  u(x,y,z,t).  The  combined 
perturbations  in  refractive  index  are  6 c/c  + u/c.  The  second  term,  which 
has  been  generally  neglected,  does  become  important  under  two  rather 
special  circumstances:  (a)  In  sound  transmission  along  a deep  downward 
loop,  and  (b)  in  the  dififieAence  between  transmission  from  a source  at 
point  A to  a receiver  at  point  B and  transmission  from  a source  at  B 
to  a receiver  at  A. 
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* I.  PROPAGATION  VELOCITY 


I ,e t r.(x,y,z,t)  designate  the  upward  particle  displacement  due 
f to  internal  wave  motion,  and 


» 


» 


u,  v,  w = 3 5 (1) 

tlie  three  components  of  internal  wave  particle  velocities.  The  perturbation 
in  sound  velocity  arising  from  the  vertical  displacement  is 

6c  = 3 c •<;  (2) 

z P 


» 


» 


I 


where  3 c 
z p 


is  the  potential  sound  velocity  gradient. 


Thus 


6c  ± u 


(3) 


is  the  i.w.  perturbation  in  propagation  velocity  in  the  ± x-direction, 
ignoring  the  small  tilt  of  sound  rays.  The  refraction  parameter  is 


» 


V 


6c  ( u 
c " c 


(4) 
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Cot'  a canonical  ocean 


l 


Lq.  {*')  can  he  written 
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1 1 . SPECTRA 


We  first  refer  to  the  previous  results  for  the  effects  of 
vertical  particle  displacements . from  MZ  (115)  and  MZ  (118)  we  have 
(lfl'>)  = (|X|2)  . Then  from  (MZ66) 


<♦*> 


wi  th 


R n 

= ~ q?nQb  f dx  nec^O  Y j-1  f dui(u j2-^)-'5  S(ui,j;z) 
) A J J 


II 

Y j del  S(u>,  j ) = > 


We  now  generalize  (2.1)  to 

<*v>  = ...  | ... 


-I 


..  S++(oi,j;z) 


+ 


and  similarly  for  (>)>  <fi  ),  (ft),  where 


L 

j 


dio 


in 


.++, 


(id,  j ) 

/ + +\  ,6c 
= (4  P > = <- 

t <“  U->  + 2 

c ' vc  c' 

j ) 

= <pV>  = < 

> + < 

> - 2 

(w,  j) 

= <mV>  = ( 

> - < 

> . 

,6c  u 


(1) 


(2) 


(3) 


(4) 


It  is  convenient  to  denote  the  contributions  to  the  spectrum  by 


s/«£ 

«£) 

s/!d 

s(—  -\ 

• 

b\  C 

c/’ 

S(c  c^’ 

\ c c / 

(5) 
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SO  tll.lt 

and  similarly  for  S and  S*  , usirp.  th»  sipn  r<,nv**nt  ion  in  (?.<*).  !n 
ao  r lance  with  the  Carrot  t and  Mirnk  internal  wave  mo<iel  we  write  (in  the 
not  it  i ■ .11  ot  MZ) 

s(— — ) = <*£«£>  gU>  h(  5 ) 

V c c / c c r,f, 

s(?  ?)  * <?  r>G™<“>  H(i'  <6> 

s(rf)  * H<1) 
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A.  The  lie  variance  for  internal  waves 

The  function  I is  to  allow  for-  the  phase  relation  between  u and 
C . We  shall  demonstrate  that  for  free  internal  waves  u and  e are 
in  fact  in  quadrature,  hence  1=0  . But  it  is  to  be  expected  that  for 
forced  (prowing)  internal  waves  this  is  not  the  case,  and  then  0 s I i 1 . 
In  the  latter  case  there  would  be  much  interest  in  measuring  the  uC 
and  ur,  = uw  covariances,  as  a possible  indication  of  the  vert1’ cal 
transfer  of  horizontal  momentum  associated  with  internal  waves. 

Constant  n.  The  simplest  demonstration  is  for  this  case;  the 
progressive  wave  solution  gives3 
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t 


t 


t 


w. 

— = sinO  cos*  cos(k  x+k  z-uit)  + sin0  sin*  sin(k  x+k  z-wt) 
u)a  x z ^ x z 


r,  = a cosO  sin(k  x+k  7.-<ot) 
x z 


where 


k x = kH  cos*  , 


k = k„  sin*  , 

y H 


k = k..  tanO 
z H 


But  there  is  evidence  that  the  waves  are  very  nearly  standing  waves  (equal 
energy  flux  up  and  down);  furthermore,  McComas4  has  demonstrated  that  if  this 
is  not  the  case  there  is  strong  nonlinear  interaction  to  make  it  so.  To 
obtain  the  standing  wave  solution  we  reverse  the  sign  of  both  0 and  kz 
and  add: 
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r 


— = -2  sinJ  cos<)>  sink  z sin(k  x+wt) 

i.),i  z x 

III  . 

a?  — sini)  sin<|>  sir*  z cos(k  x-ut  ) , 

u>  z x 

r - 2a  cos0  cosk  z sin(k  x-ut) 

s Z X 


for  progressive  waves 


C 


0 


r 


e 
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(uO  = h a2ui.  sin0  cos0  siw}> 

' ' in 

and  this  vanishes  for  a horizontally  isotropic  field,  as  assumed  GM7555  , 
since  (sin$>  = 0.  Tor  standi ng  waves  (uC)  is  proportional  to  (sinv) 
and  (cos<(>)  and  so  again  vanishes  for  a horizontally  isotropic  field. 

But  in  addition,  u and  C are  in  quadrature  vertically,  and  this  again 
leads  to  cancellation. 

Suppose  the  upward  traveling  waves  (downward  energy)  have  an 
amplitude  a,  and  the  downward  traveling  waves  ra.  Then  it  follows  that 

1 = (1-r)  (sin$) 

which  vanishes  in  the  case  of  vertical  and  horizontal  isotropy. 

General  case.  One  night  inagine  that  the  foregoing  result  is  a 

consequence  of  the  assumption  that  n is  constant,  or  of  tne  WKs 
approximation  in  which  n varies  only  slcwly.  It  is,  hcwever,  quite 
general,  let  ((k,«;z)  be  the  Fourier  component  of  vertical  displacement 
corresponding  to  an  internal  m»e  with  horizontal  waeenamber  k and 
frequency  a*.  Hub 


C 


(k,w;z)  = [d?x  j dt  e 1^X  <1’t*  r.(x,t> 


and  (the  complex  q vantity)  c(k,u>;z)  satisfies  the  equation 


L2  + k2  c(k,u,;z)  = 0 

V Z U)2-U)f  / 


for  any  buoyancy  frequency  n(z).  The  boundary  conditions  normally  used 
are  that  C should  vanish  at  the  surface  and  on  the  bottom  of  the  ocean, 
.'he  quantity  of  interest  to  us  is  (uc)  . From  the  equations  of  motion,  we 


u(k,«;z)  = k - i sz**)  3Z  c(k,w;z) 


where  e _ is  a unit  vector  pointed  vertically  upward.  It  is  then  easy 
to  show  that 


<u->  = V (A-  k » | .(k,w;z)|2> 

k z 


♦ V e xk  I*('~(k,«;z)  a *(k,w;x))) 

k*  z z 


»cw  when  averaged  over  a vertical  distance  small  caspared  to  the 
distawe  over  witui«.h  varies  appreciably.  the  first  ter®  Is 

zero-  Ant  fro®  the  differential  equation  C9)  we  can  readily  derive  the 
"conservation  law"  that 
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r 


r 


r 


I !!>('■-  ( k ,u) ; ?. ) .*  » (k  ,u;z)  ) 


is  a constant  independent  of  depth.  Since  in  particular  it  vanishes  on 

-*  . 

the  bottom,  it  vanishes  everywhere.  Thus  (u')  is  zero  for  any  buoyancy 
frequency  n(z),  without  approximation,  provided  only  that  (9)  holds. 


r 


3.  ML AN -SQUARE  PHASE 


Substitut ion  of  (?.b)  into  (2.1)  leads  to 


<4><0 


= 2n_  1 — 2 q2b(j**)  f dx  sec20  (( — ) ) I dui(u  -.  - i ' w 6 (u) 
u).^  j \ c / ] L in  r.r. 


2 f 


in 


plus  two  other  terms  involving  ((u/r)2)  and  r'*)J,  • and  (<^c/c/<  i/c)) 


and  Gu_  . Here  ( j~  * ) = T ]_l  H(j)  is  the  i.w.  weighted  average  oi  j-1 
The  w- integrations  con  be  carried  out  explicitly.^  Writing 


wf  = ut2  ♦ n2  tan20  = ug  D2  , 
L in  in  ’ 


y = w?  /u i2 
in 


we  have  for  the  three  w-integrols 


1/D 


l | dy  xfety 


2/1  P2-l  . D+l  \ 
it  ( D2  2D3  °R  D- 1 / 


'x 


or . / n' 
in 


_2_ 

3n 


dy 


1+y 


/(  i-y ) ( L-Dy) 


= A /._L 

3n  V D2 


2D3 

3D2  + 1 


, D+l  \ 
2D3  °g  D-l/ 


y (2) 


1 *>\ t ! 


Dy 


2 D3 


(1) 
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• he  x- i ntegrut i on  in  Lq.  (J.l)  i » . »r:  in  general  been  don*-  numeric  a i ly 
by  integration  along  the  ray  path.  It  can  he  done  analytically  in  two 
special  cases,  as  follows. 


A.  The  axial  approximation 

Here  L>=1,  0=0,  n n|  .We  note  that  the  second  integral 

(involving  (uu) ) becomes  logarithmically  infinite.  This  i treated  I r. 
the  Appendix.  The  result  then  is 


( 3) 


where  the  subscript  1 denotes  the  sound  axis  value,  and  similarly  for 
(<{>  ) and  (4l  + <0  . Here  R(  = 21  kn.  is  the  range  of  an  axial  loop. 

The  logarithmic  singularity  is  weak,  and  unless  we  are  within  100  m of 
an  axial  loop,  the  relative  contributions  of  the  three  terms  are  determined 

by  the  averages  ('  )))  . from  liq . (2.8),  setting  n5/n  = e-i,  the 

4 o 

relative  magnitudes  are 

24xio-8»e''*  = 1.2*10~8 
0 . 10x10-8,  e“  1 = . 'J4*10_8 
4h-2I  0 , 06x10" 8 = .02xl0’8  I 


I).  The  turning  point  approximation 


Next  we  discuss  the  apex  approximation.  Tor  steep  rays  we  note 
that  the  major  contributions  to  the  ray  integrals  in  ($%+)  come  fro.", 
the  (upper  and  lower)  apexes  of  the  ray.  Near  an  apex  located  it  <x,z), 
tiie  equation  of  the  ray  is 


z(x)  = z - (x-x)“V2 


where 


- 1 - 


= Ya  |l-cxp(z-zj )/b | 


Ya  = 1 . 14*  10-/>  km'1 


( KZ8S ) 


is  the  radius  of  curvature  of  the  ray  at  the  turning  point.  With 
this  expression,  we  may  approximately  evaluate  the  ray  integrals  and 
obtain  for  the  contribution  of  each  apex 


<*v>a_  ■ <r'>  |<(rf>  (f  f * <(?f>  * ? <(Tf)  > (r)!  <“> 


where  n denotes  the  value  of  n(z)  at  the  apex:  n = n(z)  . Thus  the 

('5(/c)‘  and  (<5c/c*u/c)  terms  get  their  major  contributions  from  upper 

n o 

turning  points,  while  the  (u/c)'  has  comparable  contributions  from  both 
( in  fac t somewhat  larger  from  the  lower  turning  point  since  the  radius 
ot  curvature  iJl  is  larger  there). 

C.  Reciprocal  transmissions 

Thus  in  any  one-way  experiment  other  than  for  a ray  with  a single 
deep  loop)  the  effect  of  particle  velocities  can  be  ignored.  If,  however. 


11 


w<‘  i mult  itneous  ly  transmit  in  opposite  directions,  then  the  three  term.-; 
•in  hr  separately  evaluated  by  constructing  the  three  combinations 

*4  j(<  4>+)  + (4>  4>  ) + 2($  $ )J 

'**  |(l  $ ) t (li  <f>  ) - 2(<{>  4>  )j  (5) 

V,  |<<fV>  4V>J 

respectively.  The  implications  of  thic  are  considered  in  Section  IV. 

D.  Worcester's  reciprocal  experiment0 

This  corresponds  to  a lower  turning  point  at  1.5  km,  and  a rather 
small  scale-depth,  B «s  0.75  km.  Thus  n/n^  = e-;!,  and  the  three  terms 
in  Cq.  3.4  are  of  a magnitude 

<<J+<{>+>  ~ jo. 44  + 0.10  + 0.42  i}  x 10-8 
Suppose  that  1=0,  then  we  have 

( i}>+  4>+  ) ~ 0.54 

( <t>  t"  ) ~ 0 . 54 

< $+  4>“  > ~ 0.34 

((<J>+-  <p~)2)  ~ 0.40 

This  is  consistent  with  some  preliminary  results  that  (4I+4|+)  , 

(<p  <P  ) and  )?)  have  comparal>le  magnitudes. 


12 
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IV.  KLC IPROCAL  CO  VARIANCE 


An  experiment  involving  rociproc.il  shooting  will  yield  not  only 
the  mean-square  quantities,  but  also  the  frequency  content.  To  interpret 
such  a spectrum  we  shall  want  to  reverse  the  integration  order,  and  do 
the  x- i ntegration  first.  Let 


< 4>  4>  ) 


[ du  E++(w) 
J 


« 


r 


f 


t 


* 


etc.  Then  the  last  term  of  the  Lq . ( 1.1)  can  be  written 


*4  I du>  LE++U)  - E ( ui ) 1 
1 


o 


dx  sec'0  (w^-w2) 


E(w)  = - ( i'1)  q;'BR  — r— 
tt  ' 1 u> . n 

i n o 


G _ ( u> ) 0(n-oj)  0(o)-oj  ) 
uC  L 


where  0(x)  is  a step  function,  0(x)  = 1 for  x > 0 and  0(x)  = 0 
for  x < 0.  These  0 functions  restrict  the  part  of  the  integral  on 
the  ray  which  can  contribute  to  the  spectrum  for  any  given  frequency.  At 
any  position  x along  the  ray,  the  only  frequencies  which  can  interact 
with  the  ray  (within  the  stationary  phase  approximation)  are  those 
sat isfy ing 


n(z(x))  > oj  > 4o?  + n?(z(x))  tan20(x) 

in 
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r 


»• 

r 


r 


t 


i tin  wlii-n  tin'  t'.iy  i:  sleep  ('•  is  large)  I requen<  i<".  near  in.^  are 

i'xc  luded . But  when  one  is  near  the  ray  apex  (n  large  and  9 small) 

a LI  frequencies  between  to. ^ and  n = n(z)  are  allowed. 

One  may,  if  one  wishes,  pick  out  different  parts  of  the  ray  by 

forming  the  difference  between  the  spectrum  at  two  nearby  frequencies. 

If  Aui  <<  u>  is  a small  frequency  difference,  then  for  frequencies  well 

above  o>.  we  have 
in 

2 2 

AClui)  E ( uit  Aui ) - E(m)  = - — ((~~\  ) (-r) r~ — * Aw/l/^1)  ; 

f ' r c /o  J n40)  v / 

o w=n 


V 

thus  Ai;(i.i)  receives  all  ot  its  contributions  from  that  part  of  the  ray 
where  n(z)x))  « u> 

# 

t 

t 


I 


f 


r. 
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V.  MOMENTUM  FLUX 


We  may  wish  to  interpret  the  results  of  reciprocal  transmission 
without  relying  so  heavily  on  internal  waves.  To  this  end  let  us  express 
the  measured  quantities  in  terms  of  the  spectrum  of  momentum  flux  in  the 
ocean  without  making  any  assumptions  regarding  the  origin  of  this  spectrum. 

Let  " designate  averages  along  a ray.  Thus ' 


G(t)  = R-1 


ds  u(s,t) 


£<t) 


R-1 


ds  £;(s,t) 


We  may  regard  u and  £ as  the  measured  quantities  for  they  are  so  closely 
related  to  the  measured  travel  times;  iet 


f 


f 


ft 


T~  = - (R/c)  G+  = - (R/c2)  [6£±G] 

designate  the  departures  from  the  mean  travel  time  (averaged  over  several 
inertial  periods).  Then 

Sc  = 3 c •£  = - (c2/R)  *5(T++T‘) 
z p 

G = - (c2/R)  ?s(T+-T') 


I 


(A<  tu.il  lv  the  vortical  gradiont  of  potential  sound  velocity  should  b" 
proper  ly  depth  weighted.) 


Conn ider  the  time- lagged  covariance 


c ( i ) = <u(  t ) r,(t-T )) 

Uf. 


Can  wo  manipulate  C to  get  some  measure  of  the  momentum  fl;  < (uw)  ? 
i'irst  wo  note  that 


At  ttii:'.  point  lot  us  simplify  the  geonw  try  to  horizontal  rays.  The 
switching  to  center  of  mass  and  relative  coordinates 


. i ’ . " ) , 


o = 


I __  t . M 


We  have 


R R-o/2  o R+o/2 

f f 

ds 

o+o/2  -R  o-o/2 


f 

1 

ds’ 

ds"  = 

do 

ds  + 

do 

. 

. 

. 

R o 

[ do( R-o ) + j do( R+o  ) 
o -R 


The  second  line  follows  if  the  integrand  is  a function  only  of  o, 
not  of  s.  Under  these  conditions 


C (o,t)  = <u(s,t)  c(s-o,t-t)) 


R 

C ,(t)  = R-2  do  ( R—  | o | ) C (o,t) 

u;  j 11  u c 

-R 


C ( 0 ) 
uw 


R"2 


R 

do( R- |o I ) C (o  ,0  ) 
11  uw 


-R 


R 

= -?R~2  do  o C (a, 0) 
j uw 

o 
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riux  spectrum.  The  more  generaJ  question  is  to  interpret  the 
covariance  of  the  "measured"  quantities: 


further 


C”(t  ) = <u(t)  r.(t-t)) 


r-pU)  = f dt  elt0T  c..(,) 

uc  I ur. 


C = (u(r.,t)  S(s-o,  t-t)) 


?o  that 


dk  dw  i(ka-o)T)  . . 

e Fuc(k’w) 


r (k,oj)  = | da  dT  e c (o,t) 


J 1 


r 


R 

R 

CO 

r-;U>)  = P'2 

ds’ 

J 

ds" 

< 

dt  e 

. 

dk 

?n 

F (k  ,w’> 

u; 


P‘?  (Is  ’ ds"  elko  r (k  ,u>) 

/n  uf. 

J J 


fO  j-R 

- R-?  l’u;(k,w)  | | do ( R+o  ) e ' k°  + ; do(R-o)  e 

-R  o 


- di 


— P r(k,u>)  ? do(R-o)  cosko 
2ti  ut 


— - ~ (1-coskR)  F (k  ,co) 


The  Reynolds  spectrum  F (k,u)  = iwF  (k,(*>),  and  so  will  be  related 
7 * uw 

to  the  quadrature  spectrum  of  the  measured  quantities. 
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HW.  Munk , J.  Acoust.  Soc.  Am.  (1474)  Eq.  24. 
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APPENDIX 


The  second  integral  in  (3.2)  diverges  as  D ♦ 1.  Set 


Ar  = D2  - 1 = (n/o . )2  tan20  . 

in 


Then  the  integral  approaches 


8 , 2 
3?  log  A 


We  can  express  A as  a function  of  ranje.  Let  source  and  receiver  be  on 
the  axis  at  x = + xj  R.  Then  with 


&z  = A cosix/R 


designating  the  distance  of  the-  ray  above  the  sound  axir  , we  find' 


where  R()  = VirBe  2 = 21  km  refer1  to  tin  axial  loop  (P.  < R for  upward 
loops).  It  follows  that 


„ 5z  3ir2  B R-R0  . ttx 

" = * IT  R -R^  S1”  T 


niO  3 it  n i B R-Rn  . ux 

a r<  = .. — — u.  sm  — 

u).  4 tt>.  R R R 

in  in  o 


9 


Sett inp,  H = 1 km,  ri)  = 1.9*10”''  sec”1,  u);)  - 7.3*10”‘J  sec”1  , 

R = 21  km,  we  have 
o 


A 


9.2  rin 

Ro 


7IX 

R 


8 . 

35  log 


4R, 


'J  . 2(  R-R  ) 
o 


Thus  even  it  we  are  within  0.1  km  of  the  convergence  distance,  R-R^  = 0.1  km, 
the  effect  of  currents  is  still  relatively  small.  So  the  logarithmic 
singularity  is  of  no  practical  interesl . 

n , “ ■’j  ^ 

As  an  aside,  the  singularity  arises  because  (w  -wf ) ► (ur-wr  ) 

n J L in 

L 

and  G ~ (ui^-io?  ) ‘ , so  that  the  integrand  of  I ~ (u>?-u>?  )”!  , and  I 

uu  in  uu  in  uu 

is  1 ogari  thmica  1 ly  infinity.  This  depends  then  on  the  fact  that  GM-spectrur. 
for  u has  a J cusp  (but,  of  course,  the  integral  of  the 

spectrum  and  hence  (u?)  are  finite).  If  the  GM-spectrum  had  a (w’-ug^) 
cusp,  with  s < h,  then  the  logarithmic  infinity  for  1^  does  not  arise. 

The  evidence  for  the  exponent  is  very  weak.  Observational  evidence 

clearly  shows  a cusp  (see  for  example  Carl  Wunsch,  1975,  JGR  80,  339-343, 

Fig.  1)  and  for  an  integrable  cusp  "...  we  require  0 < s < 1 . We 
arbitrarily  choose  the  mid-point  s = J.  ."  (GM72,  p.  249).  But 
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subsequently  Cairns  and  Williams  (Part  r r , Tip.  5)  find  the  vertical 
displacement  spectrum  (which  is  a less  critical  test  than  a u-spectrum) 
in  better  accord  with  s = 0.5  than  r = 0.1  and  s = 0.9.  The 
conclusion  is  that  the  singularity  is  probably  too  weak  to  be  observable; 
othervise,  it  would  suggest  an  acceptal  le  modification  of  the  GM-spectrurr. . 
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